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Abstract 

We study an inhomogeneous decay of an unstable D-brane in the context of 
Dirac-Born-Infeld (DBI)-type effective action. We consider tachyon and electro- 
magnetic fields with dependence of time and one spatial coordinate, and an exact 
solution is found under an exponentially decreasing tachyon potential, e~^ T ^^, 
which is valid for the description of the late time behavior of an unstable D-brane. 
Though the obtained solution contains both time and spatial dependence, the corre- 
sponding momentum density vanishes over the entire spacetime region. The solution 
is governed by two parameters. One adjusts the distribution of energy density in 
the inhomogeneous direction, and the other interpolates between the homogeneous 
rolling tachyon and static configuration. As time evolves, the energy of the unstable 
D-brane is converted into the electric flux and tachyon matter. 
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1 Introduction 



Time dependent process in string theory has been intensely studied in recent years. As- 
suming that an unstable D-brane decays homogeneously, the whole decay processes, in 
the vanishing string coupling limit g s — > 0, can be described by the marginally deformed 
boundary conformal field theory (BCFT) [T]. The main results of this time dependent 
solution, referred as rolling tachyon, indicate that, according to time evolution, the energy 
density remains constant but the pressure goes to zero asymptotically. 

On the other hand, spatial inhomogeneity has been another important issue. In partic- 
ular, much works has been studied on tachyon solitons, such as tachyon kinks [21 El El El E] 
and vortices [21 H|- These solitons are interpreted as the lower dimensional D-branes on 
the worldvolume of the original unstable system. Thus in order to see the dynamical 
formation of the lower dimensional D-branes, it is indispensable to take into account the 
spatial inhomogeneity in the decay process of an unstable system. 

The rolling of tachyon, which is inhomogeneous along one spatial direction, was con- 
sidered in BCFT (HI El llOj . The late time behavior of the resulting energy-momentum 
tensor is qualitatively different from the case of the homogeneous rolling tachyon. The 
relevant components of the energy-momentum tensor exhibit singularities at spatially pe- 
riodic locations within a finite critical time. These spatial singularities were interpreted as 
the co dimension-one D-branes [HI El El- This subject was also considered in the bound- 
ary string field theory [TT] or the DBI-type effective field theory p ] H3 \ ITo" ] IT!) ] IT? ] 

E3 E01 [^ • The Ref. [12] showed that the inhomogeneous solutions with a runaway 
tachyon potential formed caustics with multi-valued regions beyond a finite critical time, 
and proposed that in the presence of caustics, the higher derivatives of the tachyon field 
blow up. For this reason the DBI-type effective action is not reliable after the formation 
of caustics, since it was proposed as an effective action for the tachyon field in string 
theory where the higher derivatives of the tachyon field are truncated [2*2] . 

Another interesting aspect in the decay process is the dynamics at the bottom of 
the tachyon potential. This process has two kinds of decay products, which carry the 
effective degrees of freedom of the original unstable D-brane, such as energy-momentum 
and fundamental string charge [21] . They are called tachyon matter [IJ [21] and string 
fluid |25j . In the tachyon vacuum, the dynamics of the system is characterized as the two 
degrees of freedom fZE\ IT*"] 12*7] . The main purpose of this paper is to explore the formation 
of these final states in terms of the inhomogeneous tachyon and electromagnetic fields at 
late time. 

Let us consider the DBI-type tachyon effective action with gauge field interactions [22] 
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described by 



S = J <P +1 xC = -T P J d p+1 xV(T)V^X, (1.1) 

with 

X = detX lu , = det^ + d ll Td„T + F IH ,), (1.2) 
Ffa, = dfj,A v -d v A^ (/i, v = 0, 1, • • -,p), 

where A^ is an U(l) gauge field, V(T) is a runaway tachyon potential, and T p the tension 
of an unstable Dp-brane. This action (jl.lj) is expected to provide a good description of an 
unstable D-brane in the case that the tachyon field T is large, and the higher derivatives 
of T are small. As we have seen in the rolling tachyon solution pQ, the tachyon field T 
goes to infinity at late time of the D-brane decay process. Thus DBI-type effective action 
describes well the late time behaviors of the process. However, once we take into account 
the inhomogeneity of the tachyon field without gauge field interactions, as mentioned 
before, DBI-type effective action becomes inadequate in a finite critical time in describing 
the dynamics of an unstable D-brane [121 EH EJ E3 E3 EH 120] • Though we include the 
constant electromagnetic fields, the singularity we encounter seems unavoidable [TU] . 

In this paper, we suggest that the roles of the spacetime dependent electromagnetic 
fields are nontrivial on unstable D-brane system. We assume that the tachyon and elec- 
tromagnetic fields depend on time and one spatial coordinate under an exponentially 
decreasing tachyon potential, e _ ' T '/ . We find an exact solution in a frame which gives 
vanishing momentum density provided by an appropriate electromagnetic fields. The so- 
lution represents periodic profile along the spatial direction and it involves the interesting 
and inaccessible regions alternatively. In the interesting regions, the solution has no sin- 
gularities in time direction, and describes the late time behaviors of an unstable D-brane 
decay. 

In section 2, we describe the calculation of the exact solution for the tachyon and 
electromagnetic fields. In section 3, we analyze the late time behaviors of an unstable 
D-brane. Section 4 is devoted to conclusion. 

2 An Exact Inhomogeneous Solution 

Our purpose in this work is to understand the late time behaviors of the inhomogeneous 
tachyon condensation in terms of an exact solution. We take a specific frame which gives 
the vanishing momentum density over all space-time. 
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Equations of motion of the tachyon T and gauge field in the action are written 

as 

d, (7 Cf d v T) + T P V^X ^ = 0, (2.3) 

^( 7 Cf)=0, (/i, z/ = 0,1,2, •••/>), (2.4) 

where Cg 1 ' and are the symmetric and anti-symmetric parts of the cofactor C^ u of 
the matrix (X)^ v in Eq. (jl.2|) . and we define 

Conservation of energy-momentum is described by 

dfjT^ = d, ( 7 Cf ) = 0, (2.6) 
where T^ u is the energy-momentum tensor. Hamiltonian density H, is expressed as 

n = it t t + u i A i - c 

WW + W T + PiPi + {WdiTf + T*V 2 det(/i) - Ao^tf, (2.7) 

where T = <9 T, /ijj = 5^ + d{TdjT + F^-, (i,j = 1,2, •••£>), the conjugate momenta, 
IIt and II*, for the tachyon T and gauge field A{ respectively, and the conserved linear 
momentum, Pj, associated with the translation symmetry are 

rr = ^ = 7 c£, 
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From now on, let us introduce an exponentially decreasing tachyon potential, 

V(T) = V e- T / R , (2.9) 

where Vo and R are arbitrary constants. We consider an ansatz for fields which live on 
the worldvolume of the unstable Dp-brane, 

T = T(x°,x 1 ), E = F 02 (x°,x 1 ), B = F 12 (x°,x 1 ), (2.10) 

and, for simplicity, turn off all other components of the gauge fields. 
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Then the only non-vanishing linear momentum in Eq. ()2.8|) is 

Pi = -^(TT' + EB), (2.11) 

where T' = d\T. Here, we choose the zero-momentum frame due to cancelation of the 
effects of tachyon and electromagnetic fields, 

EB — —TT'. (2.12) 

Additionally the determinant X in Eq. (jl.2j) under condition ()2.12|) is factored as 

- X = (1+T' 2 + B 2 )(l-f 2 - E 2 ). (2.13) 



Conservation of energy-momentum, d fJL T^ u = 0, under the condition ()2.12j) leads to an 
observation that the energy density T 00 has only spatial dependence and T 11 time depen- 
dence, i.e., 



T oo = T^HT ^ + ^+g , (2.14) 



i _ T 2 — F 2 

T U = T 11 (t) = -T p V\l i + Tt2 + B2 , (2.15) 

where we used the notations, t = x° and x = x 1 . The Eqs. ()2.14j) and ()2.15j) are rewritten 
by 

T P V(T) = f(t)g{x), (2.16) 
l-T 2 -E 2 _f fm 2 



l+T' 2 + B 2 \g{x) 
where 

fit) = ^f^it), g(x) = y/f^{x). (2.18) 

Under the tachyon potential (J2.9|) . the equation of motion for the tachyon field (J2.3|) is 
simplified as 

f T" 1 

0. (2.19) 



1 _ f 2 - E 2 1 + T' 2 + B 2 R 
Using the Eqs. (|2.1fij) . (|2.17jl . and ()2.19j) . we arrive at the following important results, 

E = E(t), B = B(x). (2.20) 



5 



The derivation of this equation is rather technical and therefore is recorded in Appendix. 
The factorization (|2.16|) implies 



f(t) 



Rf(t) 
f(t) ' 



T'(x) 



Rg'jx) 



X21) 



From the relations (pTH?j) . fl22H, and (j2~TT]) we get 



E(t) 
B(x) 



aT 
1 



r 



a 



aRf 
ag ' 



(2.22) 



where a is an arbitrary constant. Inserting the expressions (j2.21|) and ()2.22|) into the 
Eq. (|2.17[) . we obtain the first-order differential equations for f(t) and g(x), 



+ 



R 2 (l + a 2 ) £R 2 (l + a 2 ) 
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a 2 g 2 



a 2 g 4 



0. 



0. 



(2.23) 
(2.24) 



#2(1 + a 2 ) £R 2 (l + a 2 j 
where £ is a positive constant. Solutions for the equations (|2.2Hj) - (|2.24jl are given by 



f(t) 



±\/£ sech 
±a/£ sec 



f t- c x 

\rVi 



+ a" 



a{x - c 2 ) 
RVl + a 2 



(2.25) 
(2.26) 



where c\ and c 2 are integration constants, which represent the translation symmetries 
along time and spatial directions respectively. Of course, the expressions f!2.23j) - ()2.2(i|) 
satisfy the 2-component of the gauge equation ()2.4j) . 



9 E = f 2 B>. 

Substituting the solution ()2.25|) into the Eq. ([2.5)1 . we find 



(2.27) 
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T P V 



1 + T' 2 + B 2 )(l — f 2 — E 2 ) 



Finally we obtain an exact solution for the tachyon field by inserting the expression ()2.25|) 
into (I2~TB1 . 



T(t,x) =R\n 



T P Vo 

7 



cosh 



\rVTT~c 



cos 



/ ot(x - c 2 ) 
Uvi + a 2 



(2.28) 
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This solution is characterized by two parameters, 7 and a. We will investigate the roles 
of these parameters in section 3. At first glance this result seems to be unnatural since 



cos 



/ a(x-c 2 ) 



T(t,x) has the periodic divergencies in the limit 

of cosine function. Actually in the spatially periodic regions at the initial time t 
which satisfy 



due to the property 




cos 



a(x — c 2 ) 



< 



cosh 



rVT- 



< 1, 



(2.29) 



the corresponding tachyon field T(t = 0, x) is negative. In these regions the DBI-type 
effective action does not provide a good description for the dynamics of an unstable D- 
brane as we explained in the section 1. In order to describe the late time behaviors of the 
decay process of an unstable D-brane, we restrict our interest to the spatially periodic 
regions which correspond to the large positive value of tachyon field. We will describe the 
details in the next section. 



3 Late Time Behaviors of the Decay Process 

It was observed in the previous section that there is an exact solution ()2.28|) for the 
exponentially decreasing tachyon potential in momentum zero frame. Our purpose in 
this section is to analyze the solution (|2.28J) in superstring theory. Since a total charge is 
conserved (we will mention later in detail at subsection 3.3) and the tachyon potential has 
Z 2 -symmetry under T — * — T, and runaway property V(dboo) = 0, we employ a tachyon 
potential composed of two parts, 

V(T) = V e"i = I t • (3.30) 

(V e^ T < (II) 

Tachyon profile is read from the Eq. ()2.28|) in the regions (I) and (II) by choosing the 
appropriate integration constants, C\ and C2, 

W cosh ( Jt t0 ,) cos (- 




(I) 
(II) 



where xi(xn) represents the spatial coordinate belonging to the region (I) ((H)), t is some 
large value introduced to figure out the late time behaviors of the inhomogeneous fields. 
The ranges of x\ and xu are given by 

X 4m+l — X l — X 4m.+li 

x± m+1 < x u < xt m+1 , (m = integer). (3.32) 
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where x^ m+1 and x^ m+1 denote 



x 



± 

4m+l 



V2VT 



a 



x 4m+l 



V2VT 



a* 



a 



± cos 1 
± cos -1 



7 



T P V 

7 



+ 



(4m + l)7r 



(4m + l)ir 



, (0 < 7 < W). (3.33) 



There are periodic regions, where the decay process of the unstable D-brane is not de- 
scribed well by the solution (j3.31J) . referred as inaccessible regions, 



"^4m+l ^ -^inaccessible ^ *^4?7i+l' and 2'4m+l ^ -^-inaccessible ^ -^4m+5 - 



(3.34) 



To illustrate the solution (j3.31|) graphically, we draw two figures for the tachyon po- 
tential V(T) and tachyon field T(t,x) in Fig.l. The arrows in Fig.l (b) represent growing 
tachyon field as time elapses. The corresponding tachyon fields spans the ranges at initial 
time t = 0, 



To <Tj{t,x) < T u 



-T max < Tn(t,x) < -T , 



(3.35) 



where 



T = In cosh 



to 



T 



v^ln 



v / 2v / TTo ! 

%v 



cosh 



7 



to 



V2VT 



Electromagnetic fields in Eq. ()2.22|) take the functional forms 



£(t) = aT(t) = <( v/IT ^ , V f v/ !t? / \ ,11 , (3-36) 



1 tan , =S , - f , (I) 



= __ T '( X ) = ; v; . (3.37) 

« — ' to QXn -1- si an 



The configuration for the tachyon field ()3.31)1 represents the time evolution of the spatially 
periodic profiles governed by two parameters, 7 and a. 7 adjusts the distribution of the 
energy density, while a is the scaling parameter of time and controls the period. We 
investigate the roles of the parameter a. 
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! 

■ 

(I) 



T(x,t) 
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T 




XI 
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1 
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XI 



Xi 



(II) 

I 



\T_ 



Figure 1: (a) Shape of V(T) = e ^ for \T\ > T . (b) Profiles of tachyon field T(x,t). 
Shaded regions represent the inaccessible regions. 



3.1 



a 



case : Homogeneous rolling tachyon 



When a goes to zero, the period of tachyon profile approaches infinity Thus the corre- 
sponding solution in Eq. (j3.31j) describes the homogeneous rolling tachyon, i.e., the spatial 
inhomogeneity is neglected. In the region (I) for m = in Eq. (J3.32j) . the Eqs. (jH.Hlj) . 
()3.36|) . and (|3.37j) are given by 



T{t,x) 


= V2\n 


B(x) 


(a 

= tan 1 - 

V > 


f(t) 


= tanh ^ 


T'(x) 


= E(t) = 



r p v 

7 



cosh 



V2 



cos 



axi 
V2 



t + tp 
V2 



(3.38) 



with 



— v 2 cos < axi < v 2 cos — - — . 

\ 7 / \ 7 / 

This analysis is also applicable to the case of the region (II). It is well-known that when 
the tachyon and electromagnetic fields in the DBI-type effective action (jl.lj) depend on 
only one spacetime coordinate, the results of equations of motion for the given system 
show that all the electromagnetic fields are constants EE]- In this limit, the period 
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of magnetic field B(x) reaches to infinity but its amplitude remains finite. Therefore 
the results in Eqs. (J3.38)) represent a homogeneous rolling tachyon with almost constant 
magnetic field and constant energy density T 00 = T p 2 Vq / j. 
The pressure and tachyon matter density are 



T n (t) =T 22 = - 7 sech 2 



(. 



t + to 



n T (t, x) = -Z- tanh (j^) • ( 3 - 39 ) 

In the limit of t — > oo, we obtain the pressureless matter with constant energy density 
and tachyon matter density 



3.2 a —> oo case : Static configurations 

On the other hand, when a goes to infinity, x^ m+1 and x^ m+1 reach to fixed finite values, 



x 



± 

4m+l 



X 



liu+1 



V2 
V2 



± cos 1 



1 



± COS 



-1 



r p V 

7 
T P V 



(4m + l)7r 



(4m + l)n 



(3.40) 



From the Eq. ()3.31|) we can easily check that the time dependence disappears and the 
configuration of tachyon becomes a static solution in this limit. The expressions of the 
Eqs. (13.31}) . (|3.3fjjl . and ()3.37|) at finite time t « a in the region (I) are given by 



T{t,x) 

T\x) 
E(t) 



In 



T P Vo 

7 

tan ( - 



cos 



( 



X\ 7T 

V7i~ 2 



X\ IX 

f(t) = B(x) = 0, 



(3.41) 



where 7 adjusts the distribution of energy density for the given static configuration. In 
the limit 7 — > 0, most of energy is localized at xf m+1 . The energy stored in half period of 
one cycle in the limit is obtained by 



f 1 dx+ f ° dx\ T 00 {x)\^ = 2V2T p V , 

J-xo J x x J 



(3.42) 



where Xq = \/I + a 2 7i 7 'y/2a. Since our exact solution is valid only for large T, the decent 
relation 1)3.42)1 is approximately correct. 
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As we have worked in the above subsections 3.1 and 3.2, a is the only parameter 
which governs the time scale and period along spatial direction for a given 7 in the 
solution (J3.31|) . Thus a is an interpolating parameter between the homogeneous rolling 
tachyon (a = 0) and the static configuration (a — > 00). 



3.3 Late time behaviors 



As time elapses, the tachyon profiles (|3.31|) in the regions (I) and (II) grow up along 
the arrows in Fig.l (b). In the interesting regions ()3.32j) . all physical quantities can be 
expressed explicitly, and the unstable system evolves without singularity in our system. 
Non-vanishing components of energy-momentum tensor are given by 



T 



on , 



•J' 



T (t) 



T 22 (t,x) 



, 2 / t + t 

-7sech 



1 - -i-r tanh 2 ( 



, (I) 
for (I) and (II), 

t+tp 



l+a 



tan' 



1 - tanh 2 ( Jm A + ^ tan 2 ( 



mil 



V v / 2v / l+c? 



(3.43) 

(3.44) 
(I) 

(II) ' 

(3.45) 



Since the momentum density flow is zero in the interesting regions of the solution ()3.32j) . 
the initial energy density distribution is not changed in time direction. As time goes to 
infinity, pressure along the inhomogeneous direction, x, goes to zero exponentially |24] . 
T 22 component depends on t and x-coordinates for finite a. However, in a — ► limit 
(homogeneous rolling tachyon limit), T 22 and T 11 share with the same behavior in time 
direction. 

The solution ()3.31|) also provides the expression for the electric flux density which 
satisfies Gauss constraint djl 1 = and the gauge equation (|2.27jl . 



07 



IT 



0TL7 



07 



tanh 



Wo 



tanh 



sec 



2 (~ 



etxi 



This expression denotes that absolute values of the electric flux and tachyon matter den- 
sities increase in the D-brane decay process. However, total string charge accumulated in 
the interesting regions is conserved in one period ; 



2 f axn 



(I) 
(II) 



(3.46) 




dx + 



'4m+l 



dx U 2 (t,x) = 0. 



(3.47) 
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Combining the Eqs. ()3.36|) . ( |3.43|L and (J3.46J) . we obtain the following two relations, 

( t + t 



l — T — E — sech 2 - u , (3.48) 
Ul + U 2 T = H 2 tanh 2 ( JJ^ ) , (3-49) 



where the Hamiltonian density 7i is given by 



H = JU 2 2 + II 2 , + T 2 V 2 (l + T' 2 + 5 2 ) = T 00 (x). (3.50) 

Since there is no singularity in interesting regions in time direction, the obtained 
solution describes safely the decay processes near the tachyon vacuum (V — > 0) in t — > oo 
limit. As time goes to infinity, the time dependent electric fields in the regions (I) and 
(II) become constants with opposite sign, 

E(t)\ T ^ = af{t)\ T ^ = { g+jf (I) . (3.51) 

These relations in tachyon vacuum (V — > 0) reproduce the well-known expression j2nH2E]; 

1 -T 2 - E 2 -^0. (3.52) 

This result leads to an intriguing observation. As we have discussed in subsection 13.21 the 
system gives a static configuration in the limit a — > oo. At initial time the electric field 
is equal to zero in Eq. (|3.41[) . As time evolves to infinity (t/a — > oo), the electric field 
becomes critical and T is suppressed to zero, 

\E(t)\^l, |T(*)l-0. (3.53) 

As we have seen in Eq. (|3.5U|) . the energy density is composed of three parts, such 
as string flux density (II 2 ), tachyon matter density (ITy), and tachyon potential energy. 
As time evolves, the contributions from II2 and ITt increase, while the contribution from 
the tachyon potential decreases. Finally the unstable D-brane disappears at the tachyon 
vacuum (V — >• 0). The resultant energy density in the tachyon vacuum is composed of 
two parts [2131211] 5 

n 2 = n 2 + n 2 -, (3.54) 

where 

( q 7 gec 2 ( axi A (j\ 

U 2 (t,x)\ T ^ 00 = aU T (t,x)\ T ^ 00 = { ^ \f^_ 2 {\ ■ ( 3 - 55 ) 
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4 Conclusion 



We have investigated the spatially inhomogeneous decay of an unstable D-brane in DBI- 
type effective action. We found an exact solution under an exponentially decreasing 
tachyon potential. The resulting solution involves the periodic inaccessible region along 
the inhomogeneous direction, while the behavior in time direction is well-defined. The 
solution is governed by two parameters, 7 and a. 7 adjusts the distribution of energy 
density, and a is an interpolating parameter between the homogeneous rolling tachyon 
and the static solution. 

It is well-known that the inhomogeneous rolling tachyon with a runaway type tachyon 
potential forms caustics with multi- valued regions beyond a finite critical time. After the 
critical time the unstable system may not be described by DBI-type tachyon effective 
action. However, as we have seen in section 3, it was possible to describe the late time 
behaviors of an unstable D-brane in the interesting regions due to the nontrivial roles 
of the spacetime dependent electromagnetic fields. Therefore our solution may open a 
possibility to find the caustic free tachyon field solution in a specific setting with spacetime 
dependent electromagnetic fields in tachyon effective field theory. 

As time evolves, all physical quantities are well defined and go to tachyon vacuum 
(V = 0) without developing further singularities in the interesting regions. Electric flux 
density, which is proportional to the tachyon matter density with a constant ratio a, 
increase in magnitude, but have the opposite signs in the region (I) and (II). They finally 
reach to space dependent finite configurations. As a result, the energy stored in the 
unstable D-brane at the initial stage is converted to that of the string fluid and the 
tachyon matter. Since two interesting regions in one cycle (See Fig.l (b)) go to the 
different vacua (T — > 00 and T — > — 00) in t — > 00 limit, inaccessible region between 
them contains a topological kink which seems to be interpreted as D(p — l)D(p — 1) or 
D(p- l)FlD(p- 1)F1 . 
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Appendix: Derivation of E = E(t) and B = B(x) 

In this appendix we present detailed derivation of E(t) and B(x). From the Eqs. (|2.16j) . 
(EHTJ) . and we obtain 

E 2 = l-f 2 -RT + Rf 2 p(x), 

B 2 = -(l + T' 2 + RT")+Rg 2 a(t), (A.l) 

where 

p(x) = ^, a(t) = j 2 . (A.2) 

Differentiations of the two equations in Eq. (|A.1J) with respect to x and t yield 

= i?/ 2 p'(x), 

2BB = Rg 2 a{t). (A.3) 



Product of the two equations in Eq. (JA.3j) gives 

f 3 a 3 

E' 2 = B 2 = -U-Jcr, (A.4) 

where we used the Eqs. f!2.12jl . (J2.21j) . and Bianchi identity for the gauge field strength, 

E' = B. (A.5) 

The relations in Eqs. (|A.3)1 and (|A.4|I imply that S, and B are separable functions 
of t and x. The separable properties for the electromagnetic fields, E and B, in Eq. ()A.1|) 
impose on the conditions for the functions, f(t) and g(x), 

1 - T 2 - Rf = cf 2 , 

1 + T' 2 + RT" = dg 2 , (A. 6) 

where c and d are arbitrary constants. Insertion of the relations (jA.6|) into the Eq. (jA.lJI 
produces 



£ 2 5 2 = / 2 ^ 2 ( c + Rp( x ))(-d + Ra(t)) = = T 2 T' 2 . (A.7) 



; 2 ^ 2 
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The Eq. flA.7|) is separated into two parts, 

T 2 

RT = — + df 2 , 

RT" = (T' 2 -cg 2 , (A.8) 

where ( is an arbitrary positive constant. The relations between Eqs. f|A.6|) and ()A.8|) 
give 

l + f 2 + (c + d)f 2 -l = 0, 

l + ()T ,2 -{c + d)g 2 + l = 0. (A.9) 



Differentiating the relations in Eq. (jA.9|) . we find 

f ((c + d) 

= p = RiTTo = constant ' 

T" c + d .. . 

p(x) = — = - r q + q = constant. (A. 10) 

Using the Eqs. ()A.1|) and (|A.10|) . we finally obtain 

E = E(t), B = B(x). (A.ll) 
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